
Theta Modeling Technology with Mathematical Functions and Numerical Techniques 
ispective an Underlying State Variable, Theta, wherein, an investment or derivative security 
is modeled by using a single underlying state variable, such a theta variable, insured deposits 
closed. These have a positive nominal value, an actual amount of cash value. This underlying 
state variable, theta, 0, is held to follow an independent Maikov process: d9/9 = m dt + s dz. 

This asserts that the future value of 0 depends on the known present values under a 
continuous pricing constraint. As Wiener process, dz is related to dt: Az = eVAt. 

Such theta variable depends solely on itself and time to define its expected drift and 
volatility, which it redefines through the course of its life. Thus, d0/0 = m (0,t) dt + s (0,t) dz. 

For methods drawing from standard normal distribution, the log of the change of theta 
over time and/or the log of theta at exercise should have this distribution. The theta approach 
is useful where a target variable is not die price of a traded security, but it is useful there, too. 

To creating a tradable instrument for a theta variable, assigning the function, f, as the 
price of a security dependent on 0 and time. For instance, for the insured banking's variable, 
the "deposits closed" and "deposit loss" are candidates for industry's theta. For the insured 
catastrophe dollar risk, they are the "catastrophe loss" and "net statutory underwriting loss". 
Variables are created from divers theta, such 9i, e.g. 0b and 0c, correlating respective losses. 

For instance, for ©(banking: of insured deposits closed or net deposit losses) and 0 
(cat: of insured catastrophe or net statutory underwriting losses), let f(b) and f(c) be the 
respective price of a derivative security with payoff equal to a functional mapping of 0b and 
0c into the future. Let the processes of f(b) and f(c) be defined via Ito's lemma, where 
df/f = |i dt + a dz. This stands for any f(theta). 



On a continuous time basis, the change in the price of security (f) dependent on the 
banking losses is df(b) = ub fb dt + ab fb dz; (and for Gc: df(c) =uc fc dt + oc fc dz). An 
instantaneously riskless portfolio can be created from a combination of related f(bi), such that 
( ul - r)/al = (»i2 - r)/a2 = X, where r = the present spot, risk-free interest rate at timel and 2. 

Thus, for any f, being the price of a security dependent on only 6 and time, with 

df=ufdt + afdz, there is the parameter lambda, X = (u-r)/o\ which is dependent on 
9 and time, but not on the security f, estimating the market pricing of risk of 9 by stochastics. 

The theta variable's n is the expected return from f. The expected drift, mu, equals uf. 
Sigma, of, is the volatility of f(9), and either positively (df/d9X)) or negatively correlates to 9. 
If negatively correlated, volatility = -a, and df = mu dt + (-o) f (-dz). Variance is [(0^2)^2)] 
and dz is over an independent interval, dt = (T-t). Using Ito's lemma, the parameter, mu, 
relating u and the pricing function, is set as u* f = df/dt + m 9 df/d9 + !4 s A 2 9 A 2 d 2 f/d8 2 . 

The parameter Sigma is set, ct * f = s 9 df/d9. This results in a differential structure 

dtfdt + 9 df/d9 (m - X s) + 14 s A 2 9 A 2 d 2 f/d9 2 = r* f . 

This equation can be solved by setting the drift of 8 equal to (m - A. s), and 
discounting expected payoffs at r, the present spot risk-free (usu. U.S. Treasury) interest rate. 

Thus, under risk-neutral valuation, the drift of 9 is reduced from m, to (m - X s). 

To constructing a valuation lattice, such in discrete time, introducing the notions of 

delta, 8 = e^sqrKAt)), and of mu = [2*e A (r*At)]/[ 5 + 8 A -1]. 

Hence, o = ln(8)/(sqrt At). Next, setting values for At, sigma, r(t) and 9 (where, 9 = S, 
if modeling an equity security), calculating nodes of 9 at 9(tk) = [mu A k]*[ S^^w))]*! Go]. 



Substituting theta (0) for a security (S), S isomorphic to 9, affording an underlying 
random walk of w(k(w)), such that if w=(-l,-l,l...), 8(t(3)) = [mu A 3]*[ 8 A (-1)]*[ Go]. 

In lognormal world, this relates: In 0(tn) = [n * In mu] + [w(n(w)) * In 5] + [In Go]. 
This results in the equalities: In 8 = cr*(sqrt At) and dt = T/n. Substituting and letting k = n, 
such that tn = T, forms: ln6(tk) = [n * In mu] + [w(n(w» * (a*sqrt At)]/[sqrt n] + [In Go]. 

More simply, the expected value of G, E(ln G) = In mu + In Go. 

The variance of G, Var(ln G) = (In 8) A 2. 

The volatility of 6, Vol(ln6) = (In 8)/sqrt At. 

For a pathing tree, the node value mechanic, G(tn) = [mu^]^ 8 A w(n(w))]*[ Go], using 
logarithmic transform, node mechanic, In G(tn) = [h * In mu] + [w(n(w» * (In 8)] + [In Go]. 

By the Central Limit Theorem, the term, w(n(w))/(sqrt n) exhibits strong convergence 
to the standard normal distribution, N(0,1). The term [n* In mu] shows weak convergence to 
[(r - l A o* A 2)*T], hence its robust implementation is limited to the rigors of discrete methods. 

The term [In G(T)] is distributed as [(r - >/ 2 o*2)*T] + [ N*a*(sqrt dt)] + [In Go]. 

Non-log, [G(T)] is distributed as: [Go * e^r- l A *o A 2)*T) + (N*a*(sqrt At)]. 

Valuation of a derivative security (S) based upon the state variable theta, example, the 
European call option, with realizable cashflow only at T, value today of P, with the functional 
mapping, f(6o) = max [G(T) - K, 0], where K is strike price and G is held substitutable by S. 
By weak convergence, today's value, P(G), based on G at T, derived over normal distribution: 

P(G) = [Go * N{(rT+ln(Go/K))/(o*(sqrt dt)) + (1/2 o*(sqrt dt))}] - 

[Ke'X-rT) * N{(rT+ln(Go/K))/(a*(sqrt dt)) - (1/2 a*(sqrt dt))}]. 
P = e'X-rT) * E(ra)[ G(T) - K], where E(rn) = expected value under risk-neutral conditions. 



I. A. 

For any function, f, valuing a derivative security based on theta that pays off f(T) at 
time T, the expected risk-neutral value is f = e^-rT) * E(rn)(f(T)). This requires setting the 
growth rate of the underlying theta variable in relation to [m - \*<s], rather than as m alone. 

Thus, risk-neutral valuation for today's value, P(9), P(G)=f(8), of a derivative security 
paying off f(T) at time T, is equivalent to the risk-free discount over period (0,T) of its 
expected risk-neutral future pay-out. This narrow evaluation is valid for f only over the 
continuous segment (0,T), with determinable values of F(0) and F(T). Lattices which sub- 
divide this segment, are weakened, if their At-parameters, i.e. At= (T-t) with (Kt<T, are 
modeled using analytic values from (0,T) data sets. Any methodology which relies on 
convergence to a normal distribution for its valuation, for instance, or a sampling therefrom, is 
strictly consistent only for European-style derivatives, that is, having exercise only at T, but 
not continuously throughout the segment (0 r T). Also, it assumes the security can gain or lose 
value during (0,T), with the value of the security always non-negative. The payoffs of the Gb 
and 8c securities can be European, if these stem from the single terminal condition of theta at 
T: the selected theta variables are annual aggregates, they begin each year at 9=0 and end the 
year at T, 6>0, European (0,T) events. For rigorous risk-neutral valuation, strict conformity 
can only be assigned under a European-style (0,T)segment, variable and theta-based security. 
Such a theta can substitute for a continuously traded security after adjusting for the conditions 
that 8 at T = £6i, each 9i occurring and aggregating discretely over (0,T). For continuous 
trading, full data along the annual path of such theta over (0,T) are required to be available. 

A valuation function, V, for a security or derivative dependent only on theta and time: 
V(t,0)=V(t)= Vo* ^{(r- V 2 a A 2)t + (N*a*(sqrt t)], where V(0,T) is identity of 8(0,T). 
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To constructing a swap, e.g. between the insured deposit losses and catastrophe losses, 
modeled on theta(b) and theta(c) respectively, each having valuation function, f(b) and f(c) 
respectively. The value of the swap to the payer of the deposit losses, £(b), assuming the swap 
of all year-end aggregate losses: 

V = e*(-rT) * {E(m)[f(c)(T) - f(bXT)]}, where E(rn)= risk-neutral Expectation. 
Though a swap is composed of two sides, its value, V, is a single function. Thus, V is 
a single derivative instrument modeled by the expectation of the two functions, each 
respective of its own single theta variable. Consequently, this function, V, models a security 
dependent on unrelated underlying variables, theta(i). Each theta(i) follows a stochastic 
process of form: d0/e=mi*dt + si*dzi, with mi and si the expected growth and volatility rates, 
the dzi being Weiner processes, then substituting V for f, the total loss swap, V, has the form: 

dV/V = n dt + I[si*dzi], with u being the expected return of the swap. 
Component risk of the return to the theta(i), I[si*dzi] are adjusted if theta(i) are correlated. 

Brownian motion defines the change in the value of a variable as related to the 
variable's initial value and characteristic deviation, as well as to distinct random perturbations 
resonating variance over independent intervals. It is a discrete process that approaches 
continuous form when the intervals are small and uncorrelated. Pinned simulation fixes an 
initial and terminal value for the variable, then developing the value path in between. 

An expression of theta with respect to time and to Brownian motion, the life of theta 
over (0,T) and projected Brownian motion in simulation, can be related in derivational form: 
6(t,B) = 0o * e A [(r - V 2 a*2)t + a*B]. 
This above equation values without preference to risk, obtaining risk-neutral results. 
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For geometric Brownian motion, the theta variable must be lognormal in functionality 
(i.e. its natural log values must show distribution in line with a standard normal population). 
To implementing this when modeling theta of such distribution, as respective of time only: 
9(t) = 80 * e A [(r - y 2 a A 2)t + a*sqrt(t)]. 
The weak convergence by A8(t), requires only that the natural log of the change in 
theta shows a normalised distribution and characteristic variance (not necessarily a A 2=l). 

Allowing the notation, 8(t3) = 8o*e A [(r - V 2 <s*2)t + a*B], the change in d8, measured 
at the terminal values (0,T), with t=T-0, can be derived as: d8 = 8(B)dB + [8(t) + l A 6(BB)]dt; 
its partial derivative input parameters: 
8(B) = d/dB of 8(t,B) = a*8; 
8(BB) = 0^2 * 8; and 8(t)=[r - »/ 2 <y A 2]*8. 
This computes as d8 = a*8*dB+[(r - V 2 0^2)8 + (fc o*2)&\dt, and can be reduced to: 

d8 = a*8*dB + r*8*dt 
For normal theta variables, respective only to time and theta: 

8(t) = 8o * ^[(r- V2 <s*2)t + (N*a*(sqrt t)]; and 
V(0,T) as the identity of 8(0,T), V(t,0)=V(t)= Vo * e A [(r- !4 a A 2)t + (N*a*(sqrt t)]. 
This requires only that the natural log of the change in theta, hence, in V, has a characteristic, 
normal distribution. N represents a sampling off the standard normal distribution, e.g. N=8=4>. 

Monte Carlo simulation is a discrete, methodology that is based on the Law of Large 
Numbers (e.g. in large numbers of sampling sequences). The life of the security is subdivided 
into n intervals, each of length At Using ssvolatility, and m= risk-neutral growth rate, of 8: 
A8 = m*8*At + s*8*e(VAt), where each simulation run has n drawings, one per At.. 



For a multiple state 0: A9i = mi*9i*At + si*9i ei(VAt), with 9i: (1 <> i £ n). If the 9i 

i 

are correlated, implement correlation between the 9i, pik, and also between the si and pik. 
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The Figure 91 diagrams the mathematical programming functions for theta variable. 
The Figure 92 diagrams the mathematical programming functions for security based on theta. 

For further embodiment of theta, and as respective reinsurance and actuarial sciences, 
see inventor's published work, "Forecasting Expectations of Insured Depository Default and 
Catastrophic Losses". The publication further elaborates the purposes and mechanics of theta. 

Regarding the use of a put-option model for deposit guarantees, i.e. insurance against 
deposit losses, a deposit insurance guarantee is not isomorphic to an European put option, as 
any recovery of underwritten (collateralizing) assets, occurs later than the deposit loss occurs. 
Asset recovery or loss development does not occur simultaneously to default or catastrophe. 

An OAS/martingale lattice, featuring the development of recovery or loss over time, 
implemented in diagram schematic, Figure 93, an Option Adjusted Spread style lattice. 
Having an initial node, at t = 0, at which node having value at given t, expressed as e(0) = 1. 
From initial node, binary pathing, to the upper node at time t = 1, processing 1 - M(0)h(l), 
wherein M is martingale for the given h, and wherein h is distinct time interval, value at upper 
node of e(l) = 1, a state wherein no default or loss occurs; to the lower node at time t = 1, 
processing M(0)h(l), value at lower node of e(l) = E(L), a state where default or loss occurs, 
wherein E(L) is the risk-neutral expectation of L after asset recovery or loss development. 

From the upper node at time t = 1, binary pathing: to an upper node at time t = 2, 
processing 1 - M(l)h(2), value at upper node of e(2) = 1, wherein no default or loss occurs; to 
a lower node at time t = 2, processing M(l)h(2), value at lower node of e(2) = E(L), a state 
where default or loss occurs. If lower node occurs, follow path relation of lower node at t = 1. 



